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Recent work lias shown that a temperature dependence of the shear viscosity to entropy ratio, 
r]/s, influences the collective flow pattern in heavy-ion collisions in characteristic ways that can be 
measured by studying hadron transverse momentum spectra and their anisotropies. Here we point 
out that it also affects the pair momentum dependence of the Hanbury-Brown-Twiss (HBT) radii 
(the source size parameters extracted from two-particle intensity interferometry) and the variance of 
their event-by-event fluctuations. This observation establishes interferometric signatures as useful 
observables to complement the constraining power of single-particle spectra on the temperature 
dependence of r//s. 
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I. INTRODUCTION 

Recently, a great deal of theoretical attention has 
been devoted to the accurate determination of the sliear- 
viscosity-to-entropy-density ratio (??/s) of the quark- 
gluon plasma (QGP) and the hadron resonance gas 
(HRG), and how its temperature dependence around 
the transition point between these two phases might be 
probed in heavy-ion collisions. In particular, a recent 
study [1] suggested the possibility of discriminating ex¬ 
perimentally between different parametrizations of the 
temperature dependence of p/s using systematic anal¬ 
yses of the anisotropic flows and iq as functions of 
transverse momentum pt and pseudorapidity rj. In this 
paper, we present a complementary analysis of the HBT 
radii for Au+Au collisions at RHIC energies, to argue 
that the HBT radii can be used profitably to help con¬ 
strain viable parametrizations of (77/s) ( T ). 

The HBT radii are physical observables which are 
well defined on an event-by-event basis. As parameters 
that describe the size and shape of the particle emitting 
source, they fluctuate from event to event. Experimen¬ 
tal data so far [2-7] determine only the average source 
radii characterizing a large ensemble of events [8]. Else¬ 
where [9] we discuss new ideas of how to complement such 
measurements of the ensemble-averaged HBT radii by a 
determination of their variances, i.e., the width parame¬ 
ters characterizing their event-by-event fluctuations. In 
that work [9] it is also shown that the experimentally ac¬ 
cessible observables characterizing the mean source radii 
and their variances clearly track the algebraic mean and 
variance of the source radii that are calculated theoreti¬ 
cally from the individual emission functions of a large set 
of dynamically simulated heavy-ion collisions with fluc¬ 
tuating initial conditions. We here show the transverse 
pair momentum (Kt) dependencies of these mean radii 
and their variances for the same four parametrizations of 
the specific shear viscosity 77/s studied in [1]. We also 
compare the mean radii to those one would obtain from 
a single average source by dynamically evolving a single 
initial profile, constructed as the ensemble average of the 


fluctuating initial conditions of the full ensemble. 

To save time and effort, we employ for this pioneer¬ 
ing study the shortcut of calculating the HBT radii from 
the source variances of the emission function [10-13] for 
directly emitted pions only, instead of performing a 3-D 
Gaussian fit to the correlation function including all res¬ 
onance decays [13-16]. This approximation, which holds 
exactly only for Gaussian sources [12], is known to be 
sufficient to discern qualitative features of the HBT radii 
and their Kt dependences, although it is not accurate 
enough for quantitative comparisons with experimental 
data [17]. 


II. ANALYSIS 
A. Initial conditions 

In this paper, we compute the HBT radii, as outlined 
in the previous section, for central (0-10%) Au+Au col¬ 
lisions at 200 A GeV, based on charged pion correlations. 
We use the iEBE-VISHNU package [18] in modeling the 
evolution of heavy-ion collisions. For the (hydrodynami- 
cal) initial conditions of the fireball, we employ the MC- 
Glauber model [19] with p+p multiplicity fluctuations 
[18] to generate the initial entropy density distribution 
in the transverse plane. The subsequent dynamical evo¬ 
lution is performed on an event-by-event basis, meaning 
that each set of initial conditions that we generate ini¬ 
tializes one of the N ev hydrodynamically evolved events 
in our ensemble. In this paper, we take N ev = 1000. 


B. Hydrodynamic evolution and ( r]/s)(T) 

We perform the hydrodynamical evolution of the ini¬ 
tial conditions with the VISH2+1 code [20], using the 
s95p-PCE165-vO equation of state [21]. We choose this 
evolution to begin at an initial proper time To = 0.6 fm/c, 
where t = 0 corresponds to the instant at which the two 
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nuclei collide. Also, although necessary for quantitatively 
precise comparisons with experimental data [22], we do 
not incorporate here any initial transverse flow into the 
hydrodynamical evolution. 

For our study, we use for ( r)/s)(T ) one of the four 
parametrizations given in [1]: 

• LH-LQ, in which (ry/s)(T) = 0.08 for all T; 


{r]/s)(T) qgp = -0.289 + 0.288 ^ + 0.0818 ( ^ 
for T > T tr ; 


(t7/s)(T) H rg = 0.681 — 0.0594 — 0.544 ( A_ 
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and (t]/s)(T)q G p for T > T tr . 


Here, T tr = 180 MeV [1] represents the transition tem¬ 
perature between the QGP and HRG phases. The four 
parametrizations are shown in Fig. 1. 
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same final thermal 7r + multiplicity. 1 Normalizing all of 
the averaged profiles to the ideal fluid case (which thus 
corresponds to Kg t = 1) we find 

• LH-LQ: Kfit = 0.889 

• LH-HQ: K fit = 0.596 

• HH-LQ: K fit = 0.878 

• HH-HQ: K fit = 0.589 

As expected, parametrizations with a large plasma vis¬ 
cosity require a significant suppression of the initial en¬ 
tropy density profile in order to compensate for viscous 
heating and to reproduce the measured final multiplici¬ 
ties in Au+Au collisions at RHIC energies. Finally, we 
terminate the hydrodynamical evolution on a decoupling 
surface (the ”freeze-out surface”) of constant tempera¬ 
ture Xdec = 120 MeV. 


C. Pion emission function at freeze-out 

On the freeze-out surface for each event, the pion emis¬ 
sion function is defined by the Cooper-Frye integral [23- 
25] 


S(x,p) = J^P'd 3 cr{y)5 4 (x-y)f{y,p), (3) 

f(x,p) = f 0 (x,p) + Sf(x,p) 

= eipu-^/T^ + 2T2 P (e+P) /o(1+/o) - (4) 

Here, Sf is the first-order viscous corrections to the ideal 
distribution function f 0 [26, 27] for which we assume a 
quadratic dependence on p. tt^^x) and u^(x) are the 
shear stress and flow velocity profiles along the freeze-out 
surface E, respectively, p, T, e, and V are the chemical 
potential, decoupling temperature, energy density, and 
pressure, respectively, which are all constant on E by 
construction. d 3 a fl (x) is the outward pointing normal 
vector on E at point x. 


FIG. 1: The different parametrizations of ( y/s)(T ) used in 
this work, taken from [1], 

As also pointed out in [1], for identical initial condi¬ 
tions viscous entropy production generates different final 
multiplicities for different choices of ( 77 /s) (T). This ef¬ 
fect is especially severe for the parametrizations LH-HQ 
and HH-HQ where at early times, when the longitudinal- 
transverse velocity shear is biggest, ( 77 /s) (T) is very 
large. To compensate for this effect, we rescale our initial 
entropy profiles by a constant factor Kg t , which depends 
on the chosen parametrization of ( 77 /s) (T). For a given 
choice of ( 77 /s) (T), Kg t is chosen so that an ensemble- 
averaged initial entropy profile always reproduces the 


D. Ensemble averaging 

Theoretically, the HBT radii may be defined for either 
a single event or an ensemble containing many events. 
Due to statistical limitations arising from the finite num¬ 
ber of particles emitted in a single event, precise mea¬ 
surements of the full set of HBT radii and their Kt de- 


1 Since we always decouble at the same constant freeze-out tem¬ 

perature, —is a fixed fraction of the total multiplicity '//) . 

Therefore, this normalization ensures that, on average, all events 

have the same final charged multiplicity. 
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pendences for a single event are not possible [8]. Three- 
dimensional experimental HBT analyses, therefore, ex¬ 
clusively report measurements based on large collections 
of events, rather than measurements of event-by-event 
HBT radii. In order to obtain an apples-to-apples com¬ 
parison with experimental data, theoretical HBT calcu¬ 
lations must therefore, at some level, average over all 
events in the ensemble. 

This procedure, known as ensemble averaging, can be 
performed in several ways, two of which we consider here. 
According to one prescription, the initial conditions for 
all events in the ensemble are averaged before their hy- 
drodynamical evolution. We refer to this method of en¬ 
semble averaging as ’’single-shot hydrodynamics” (SSH) 
[30]. For a sufficiently large ensemble, averaging over 
many, individually fluctuating and bumpy initial condi¬ 
tions results in a smooth initial transverse entropy density 
profile, which in turn may be evolved hydrodynamically 
as a single, averaged event. This prescription eliminates 
any sensitivity to event-by-event fluctuations in the ini¬ 
tial state. We will here denote the HBT radii derived 
from single-shot hydrodynamics by R 2 j ( i,j = o,s,l). 

An alternative prescription for computing the 
ensemble-averaged HBT radii evolves each fluctuating 
event independently and averages over the entire ensem¬ 
ble only after the HBT radii have been computed for each 
event from its individual hydrodynamic emission func¬ 
tion S(x,K). This leads to the following definition of 
ensemble-averaged HBT radii: 

ev fe=x 

We denote by ( R 2 ,-), the HBT radii of the kih event, and 
we label the mean HBT radii computed according to this 
’’direct ensemble average” (DEA) by (i??-). 


E. HBT calculations 

The basic formalism describing our calculation of the 
HBT radii was presented in [8]. Once the emission func¬ 
tion (3) for a particular event has been obtained from 
the hydrodynamic output using Eq. (4), the HBT radii 
corresponding to that emission function may be defined 
by [10, 12] 


= {{xi - Pit){xj - Pjt)) s , (6) 


where = x y .-(x yL ) s , (3 = K/E k , E k = y ml+K 2 , and 

the I\ dependence of Rf :) arises from the K dependence 
of the emission function in the definition of the source 
average: 


, . _ f d 4 x f(x) S{x,K) 

' X s f d 4 x S(x, K) 


( 7 ) 


K = is the average pair momentum in the 

lab frame, K = | (p-\ +p -2 ). In this paper, we consider only 
pairs of identical particles at mid-rapidity (Y = 0). 

The dependence of the HBT radii on the pair emission 
angle <&k also allows them to be expanded in a Fourier 
series (we suppress the additional dependence of all HBT 
parameters on Kt and Y), 

OO 

Rl^ K ) = i4 0 + 2^ cos(n($K—4/„)) 

n—1 


R l\n , (8) 


where the flow-plane angle is defined by the complex 
phase of the nth-order anisotropic flow coefficient v n [31]: 

V <n#B = «e in *» (9) 

/ 0 °° d,K T K T f\ d$ K e in4>K f d 4 x S(x, I<) 
f 0 °° dK T K T JZr d^K f d 4 x S(x, K) 


In this work, we average R 2 , R 2 0 , and R 2 over § K , but we 
report the 2nd <f>x-harmonic R 2 s 2 of the out-side cross¬ 
term R 2 s (&k), since the mean value of its d^-average 
(R-os o ) vanishes for symmetric nucleus-nucleus collisions 
at midrapidity [32]. For single events with fluctuating 
initial conditions, R 2 S 2 receives non-zero contributions 

from both R 2< f\ an d on an event-by-event basis, but 
the former contribution averages to zero for a sufficiently 
large ensemble. We therefore consider only the mean and 
variance of the event-wise fluctuations of R^\ — R os 2 
for the remainder of this paper. For notational simplic¬ 
ity, we will drop the ’0’ subscript from the azimuthally 
averaged radii, and write simply R 2 :j ( i,j = o,s,l ). 

At Y = 0, Eq. (6) gives 

R 2 s = (x 2 ) s , R 2 = (x 2 ) s . (10) 

R 2 and R 2 thus depend only on the geometric aspects of 
the emission function. By contrast, both R 2 and R‘ 2 t in¬ 
volve a combination of spatial and temporal source vari¬ 
ances: 


R l = ( S o> s - 2/3 T (i 0 f ) s + ^(t 2 ) s , (11) 

R 2 os = {x 0 x s ) s - /3 T (xsi) s (12) 

In [9] we describe a method for experimentally measur¬ 
ing, in addition to the mean, also the variance of an event- 
by-event distribution of the HBT radii. As a prediction 
for such measurements, we here also compute the vari¬ 
ances of the HBT radii for our ensemble of N ev = 1000 
events, defined by 2 

1 JVev 2 

4 = jriE((4)r(4>) ■ ( 13 ) 

ev k =1 


2 The factor N 1 _ 1 makes o - ?. an unbiased estimator of the true 
variance of the event-by-event distribution of HBT radii [33]. 
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FIG. 2: Rfj vs. Kt from single-shot hydrodynamics, for ideal 
and viscous fluids with the temperature dependent shear vis¬ 
cosities (rj/s) ( T ) shown in Fig. 1. Note the almost perfect 
overlap of the green (dashed) and blue (dotted) curves in pan¬ 
els (b,c). 


In this paper, we will use and plot aij/ ( R?j ) as a measure 
for the widths of the HBT radii distributions. For the 
quantity R 2 S 2 , we define the corresponding variance to 
be 


'os, 2 


N ew 2 

■E((^.,2) fc -<^.,2>) ■ ( 14 ) 


fc=l 


III. RESULTS 

A. Mean HBT radii and their Kt dependence 


In this subsection, we present the mean HBT radii 
computed according to the two averaging prescriptions 
described above (single-shot hydrodynamics (SSH) and 
the direct ensemble average (DEA)), and study their sen¬ 
sitivity to the temperature dependence of the specific 
shear viscosity ( 77 /s) (T). 

First we consider in Fig. 2 the radii A/ extracted from 
single-shot hydrodynamics which averages over event-by- 
event fluctuations in the initial state. Panels (a) and 
(d) show that the sideward radius R 2 S and the out-side 
cross term R 2 S 2 exhibit sensitivity to the temperature 
dependence of ( 77 /s) (T) both below and above T tr . How¬ 
ever, ( 77 /s) (T) affects hadronic observables in two dis¬ 
tinct ways. First, there is the cumulative dynamical 
effect of shear viscosity on the development of flow in 
the fireball; flow probes the entire temperature history 
of ( 77 /s) (T) between the initial and decoupling temper¬ 
atures. Second, the emission function is affected by 
the viscous correction Sf to the distribution function at 
freeze-out, Eq. (4), which is controlled by the value of 


FIG. 3: Same as Fig. 2, but leaving out in Eq. (4) the viscous 
correction Sf at freeze-out. Comparison of Figs. 2a,d and 
3a,d shows that the sensitivity of Rs(Kt) and Ro S , 2 (Ar)to 
the temperature dependence of (rj/s) ( T ) below Tt r is almost 
entirely due to the theoretically poorly controlled viscous cor¬ 
rection Sf at freeze-out, and that any sensitivity to the tem¬ 
perature dependence of (77 /s) (T) during the dynamical evo¬ 
lution, i.e. before freeze-out, is strongly weighted at high tem¬ 
peratures, i.e. early times. The sensitivity to Sf at freeze-out 
of the HBT radii probing the outward and longitudinal di¬ 
mensions of the emission function (panels (b,c)) is negligible. 


( rj/s)(T ) on the freeze-out surface. This gives a con¬ 
tribution to observables such as the HBT radii that de¬ 
pends exclusively on the behavior of ( 77 /s) (T) at Xd ec 
and is basically independent of its prior history (except 
for the accumulated effect of (77 /s) (T) on the flow pat¬ 
tern whose associated shear tensor also affects the shear 
stress at freeze-out). To separate the two effects, we plot 
in Fig. 3 the HBT radii computed without the Sf cor¬ 
rection, so that only the cumulative dynamical effects of 
( 77 /s) (T) remain. Figures 3a,d show that deleting Sf re¬ 
moves all sensitivity of the sideward and out-side HBT 
radii on the shear viscosity in the hadronic phase, leav¬ 
ing only a weak sensitivity on the behavior of ( 77 /s) (T) 
in the high-temperature QGP phase. 

On the other hand, comparison of panels (b,c) in 
Figs. 2 and 3 reveals that the outward and longitudi¬ 
nal HBT radii are almost completely unaffected by the 
behavior of ( 77 /s) (T) below T tr and, in particular, receive 
no significant contribution from the viscous correction Sf 
at freeze-out. This finding is at variance with the results 
reported in Teaney’s pioneering analysis [26] of /-effects 
on the HBT radii. We note, however, that Teaney’s anal¬ 
ysis was based on a blast-wave parametrization of the 
hydrodynamic flow profile at freeze-out, rather than the 
profile from a genuine hydrodynamic simulation of the 
flow. It was noted before in Appendix E of Ref. [20] 
that Sf is very sensitive to the details of the velocity 
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K t (GeV) K t (GeV) 


FIG. 4: The ensemble averaged (R?j ) as a function of pair 
momentum Kt, for the same choices of (p/s) (T) as in Figs. 2 
and 3. This figure should be compared with Fig. 2. 


shear tensor at freeze-out, and even its sign is fragile. 
Furthermore, there are well-known conceptual uncertain¬ 
ties about the correct form of Sf, especially its pt de¬ 
pendence [27]. These uncertainties render suspect any 
observable that strongly depends on Sf at freeze-out. 3 
Therefore, we caution the reader not to trust the appar¬ 
ent sensitivity in Fig. 2a,d of R 2 (Kt) and R 2 S 2 (Kt) to 
the behavior of (p/s) (T ) in the hadronic phase. On the 
other hand, the sensitivity of R 2 (Kt) and R/(Kt) ex¬ 
clusively to the temperature dependence of (p/s) (T) in 
the QGP phase above T tI appears to be robust and un¬ 
affected by freeze-out uncertainties related to Sf. These 
high-temperature effects of (p/s) (T) on R 2 and Rf are 
not huge, but they can reach 20% (for the squared HBT 
radii) for a range of Kt values below 1.5 GeV. 

Studies such as [15, 16] have explored the sensitivity 
of the HBT radii to shear viscosity using constant (T- 
independent) p/s and varying that constant. Different 
from our approach here, when changing p/s they also 
changed other hydrodynamic parameters to ensure that 
not only the normalization but also the slope of the px~ 
spectra was held fixed. Additionally, both of these stud¬ 
ies extract the HBT radii from a Gaussian fit to the full 
three-dimensional correlation function instead of using 


3 We comment that <5/-related uncertainties may also affect some 
of the conclusions in Ref. [1] about the dependence of charged 
hadron elliptic and quadrangular flow, v 2 (Pt) and v/' (pt), on 
the behavior of (p/s) (T) in the hadronic phase. For example, we 
found that, when calculated without the 5f correction, almost 
all sensitivity of v 2 (p r r ) of directly emitted (“thermal”) pions to 
the temperature dependence of (p/s) (T) in the hadronic phase 
disappears, leaving only a (significantly weaker, but robust) de¬ 
pendence of this observable to the temperature dependence of 
(p/s) (T) in the QGP phase above Tt r . 



K t (GeV) K t (GeV) 


FIG. 5: Same as Fig. 4, but without the Sf correction. This 
figure should be compared with Fig. 3. 


the short-cut through the source variances. These differ¬ 
ences make a direct comparison of our work with their 
results difficult. Our findings are qualitatively supported 
by the work [28] (whose authors, like us, used the source 
variances method) where it was found that an increase 
in the Knudsen number (or, equivalently, p/s) resulted 
in a decrease in R 2 . The purpose of the present work is 
not a realistic comparison with experimental data which, 
as done in [15, 16], would require a simultaneous tuning 
of several additional parameters (see also Ref. [29]). Our 
goal is to check systematically the sensitivity of the HBT 
radii to a possible T-dependence of p/s , in particular in 
the context of event-by-event source fluctuations. It is to 
these fluctuations that we turn our attention next. 

We now investigate how the mean HBT radii are af¬ 
fected by event-by-event fluctuations in the initial condi¬ 
tions. To do so we evolve 1000 central (0-10% centrality) 
Au+Au events at y fs = 200 A GeV with fluctuating initial 
conditions and average their correspondingly fluctuating 
HBT radii after freeze-out, as described in Sec. IID. The 
results are presented in Fig. 4 (including the full distri¬ 
bution function (4) at freeze-out) and Fig. 5 (without the 
viscous Sf correction in (4) at freeze-out). Compared to 
Figs. 2 and 3, we see no qualitative differences. Again, 
robust (although not very strong) sensitivities to the tem¬ 
perature dependence of (p/s) (T) above T tI are seen in all 
radii, but for R 2 these are buried by a stronger sensitiv¬ 
ity to the viscous Sf correction at freeze-out (which, we 
reiterate, is theoretically not well controlled). As was the 
case for the R 2 j from single-shot hydrodynamics, the Sf 
corrections to the (i??-) are negligible when i,j = o, l. 

In Fig. 6 we plot the fractional change 
((R 2 j)—R 2 j) /B ,2 j of the mean squared HBT radii arising 
from event-by-event fluctuations, as a function of pair 
momentum. Generically, event-by-event fluctuations are 
seen to boost the T^-averaged means of the fluctuating 
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FIG. 6: The fractional change ((7?/,)—A/,) /B%j of the 
squared HBT radii due to event-by-event fluctuations, as a 
function of pair momentum Kt■ See text for discussion. 


radii by a few percent above the corresponding radii 
obtained from single-shot hydrodynamic evolution of 
a smooth averaged initial condition. Exceptions to 
this rule are the sideward radius for the HH-LQ and 
HH-HQ parametrizations (i.e., for large shear viscosities 
at freeze-out) in the range K T > 0.5 GeV, driven by 
effects from the 6f correction at freeze-out, and the lon¬ 
gitudinal radius at small pair momentum Kt < 0.2 GeV. 
The systematic effects of event-by-event fluctuations 
on the mean $ A '-oscillation amplitude ( R 2 os 2 ) are less 
unambiguous, with fluctuations increasing or decreasing 
(R 2 S 2 ) relative to R 2 s2 hr different R^-ranges depend¬ 
ing on the specific temperature dependence selected 
for (77/s) (T); the differences can become substantial 
at small Kt where the oscillation amplitude is small, 
with (i?o s 2 ) even changing sign relative to R 2 S 2 in the 
LH-HQ and HH-HQ parametrizations. Generically the 
differences between the mean fluctuating HBT radii and 
their single-shot hydrodynamic analogs are largest for 
ideal fluid evolution and somewhat smaller for viscous 
evolution. This is consistent with the idea that viscosity 
attenuates fluctuation-driven effects. 

For a better understanding of how different tempera¬ 
ture dependences of ( 77 /s) (T) affect the magnitudes and 
Kt dependences of the radius parameters and ( Rfj ) 
we can look directly at the emission functions. These are 
shown in Fig. 7a for a randomly selected single fluctuat¬ 
ing event, and in Fig. 7b for the single-shot hydrodynamic 
event with smooth ensemble-averaged initial conditions. 
Having seen in Figs. 2-5 that at RHIC energies the HBT 
radii discriminate mostly between ( 77/s ) (T) parametriza¬ 
tions that differ in the value of 77 /s at the earliest times 
and highest temperatures, we focus in Fig. 7b on the 
HH-LQ and HH-HQ parametrizations (a comparison of 
LH-LQ with LH-HQ would lead to similar conclusions). 


As is well known (see, e.g., Ref. [34]), the average ra¬ 
dial flow velocity increases approximately linearly with 
r as one moves from the top of the freeze-out surface 
in Fig. 7b (at r = 0) to its vertical part (where r = r max 
takes its maximum value). For the cases studied in Fig. 7, 
this is shown in Fig. 8 . This observation explains why, 
for increasing I\t , the region of maximal emissivity (col¬ 
ored red) moves from around r = 0 for small momentum 
pairs to r max for large momentum pairs. Fig. 7a illus¬ 
trates how this phenomenon manifests itself in an indi¬ 
vidual fluctuating event: While the regions of highest 
emissivity roughly follow the average tendencies seen in 
Fig. 7b, these tendencies are strongly modulated by fluc¬ 
tuations in the shape of the freeze-out surface and of the 
radial flow along this surface, giving rise to peaks and 
valleys of emissivity as one moves around the freeze-out 
surface at constant values of r or r. Fig. 7a also illus¬ 
trates that, in an individual event, the regions of maximal 
emissivity have a strong dependence on the azimuthal di¬ 
rection of pair emission, and that therefore, for an 
individual fluctuating event, one should expect strong 
4>#-dependence of the HBT radius parameters. 

Comparing the freeze-out surfaces and regions of max¬ 
imal emissivity shown in Fig. 7b for the two ensemble- 
averaged sources that were evolved with shear viscosities 
( 77 /s) (T) differing only in their behavior at high temper¬ 
atures, we can make several observations: 

(1) The larger specific shear viscosity of the HH-HQ 
parametrization drives stronger radial flow which causes 
the fireball to expand to larger radii and complete 
its freeze-out at earlier times than for the HH-LQ 
parametrization. This agrees with similar observations 
made in Ref. [20]. 

(2) For longitudinally boost-invariant expansion with 
longitudinal velocity component v z = z/t , emissivity re¬ 
gions centered at later times see a smaller longitudi¬ 
nal flow gradient and thus a larger longitudinal region 
of homogeneity, reflected in a larger value for Rf [12], 
Fig. 7b shows that maximal emissivity for small and 
intermediate Kt pairs is shifted to later times for the 
HH-LQ parametrization than for HH-HQ, while for large 
Kt = 2 GeV the pairs are emitted at roughly the same 
times for both parametrization. This explains why in 
Figs. 2c (and also in Fig. 4c) the longitudinal HBT radii 
are larger for HH-LQ than for HH-HQ at small and in¬ 
termediate Kt but nearly identical at large Kt- 

(3) For small Kt, the regions of highest emissivity are 
geometrically smaller in the outward and sideward 4 di¬ 
rections for the HH-LQ case than for HH-HQ, but the 
opposite is true at intermediate and high pair momen¬ 
tum. This explains the behavior of the sideward radius 
R 2 S for these two parametrizations shown in Figs. 2a. 


4 To see the truth of this second statement the reader should in 
her mind visualize the emission surface rotated in </> direction 
and look at its horizontal width when frontally viewed. 
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FIG. 7: (a) The emission function S(x,K), for three values of Kt as indicated, F' = 0 and integrated over and r/ s , for 
a single fluctuating event from the 0-10% centrality range, evolved with (?7/s)hh-lq. Colors code the emission intensity, 
normalized to its maximum value on the freeze-out surface, (b) The emission function S(x,K) for the smooth average single¬ 
shot hydrodynamic event corresponding to 0-10% centrality, for the same three values of Kt, Y = 0 and integrated over 
and rj a , using the same color code as in (a). The width of the bands arises from plotting the (r, r) freeze-out contours for all 
spatial angles (j> a on top of each other. The left and right columns of panels show results for evolution with (?;/s)hh-lq and 
(??/s)hh—hq, respectively. See text for discussion. 


The most important aspect of the freeze-out surfaces 
shown in Figs. 7b is the earlier freeze-out associated 
with larger shear viscosities at early times (high temper¬ 
atures). It reduces the emission duration and increases 
the longitudinal flow gradient at freeze-out, which re¬ 
duces the (x 2 ) s and (t 2 ) s contributions to Rf and R 2 
in Eqs. (10), (11) and is the root cause for the robustly 
smaller R 2 and R 2 values observed in Figs. 2b-c and 4b-c 
for the HH-HQ parametrization of ( 77 /s) ( T ) compared to 
the HH-LQ case. Its effect on the cross-terms variances 


(x 0 Xs) s , (x 0 P} s and (and thus the behavior of 

the cross-term radius R 2 S ) is less intuitively obvious. 

Last but not least, it is well established that the rate 
at which R 2 decreases with increasing pair momentum 
Kt is a measure of the radial velocity gradient along 
the freeze-out surface [12], Comparison of Figs. 8b,e,h 
(middle column) to Figs. 8c,f,i (right column) shows that 
evolution with (?7/s)hh-hq leads to a larger radial flow 
gradient along the freeze-out surface than evolution with 
(t?/s)hh-lq (which has smaller shear viscosity at early 














r (fm) r (fm) r (fm) 


FIG. 8: For each cell on the freeze-out surface the radial velocity is plotted against its radial position, color coded as in Fig. 7 
for its emission intensity normalized by its maximal value on the surface. Panels a,d,g in the left column are for the single 
fluctuating event, evolved with (7//s)hh-lq, shown in Fig. 7a, for the same three Kt values as indicated. Panels b,e,h in the 
middle column and panels c,f,i in the right column are for the ensemble-averaged event shown in Fig. 7b, evolved with single¬ 
shot hydrodynamics using (t?/s)hh-lq (middle column) and (?7/s)hh-hq (right column), respectively. See text for discussion. 


times). This correlates with the steeper Kt dependence 
for R 2 S in Figs. 2a and 4a for the HH-HQ parametrization 
(green dashed curves) than for the HH-LQ one (black 
dash-dotted curves). 


B. Variances of HBT radii and their Kt 
dependence 

Finally, we consider the variances of the event-by-event 
distributions of the HBT radii. We present them in the 
form of relative widths, normalized to the directly en¬ 
semble averaged squared radii, in Figures 9 (full freeze- 
out distribution function) and 10 (without the viscous 
Sf correction on the freeze-out surface). For the $k~ 
averaged radius parameters the relative widths are seen 
to be of order 10 — 15%, with little dependence on either 
Kt or the specific shear viscosity used in the hydrody¬ 


namic evolution. The event-by-event distribution of the 
out-side oscillation amplitude R 2 s2 shows a somewhat 
smaller relative width, especially at small Kt where the 
mean amplitude is small. At larger Kt, shear viscos¬ 
ity appears to reduce the relative width of the event-by¬ 
event distribution of this oscillation amplitude compared 
to ideal fluid evolution. 

Focusing on the $,R-averaged radii (panels (a-c) in 
Figs. 9 and 10) we observe that, just as was the case for 
their ensemble averages, all sensitivity of their variances 
to the temperature dependence of rj/s in the hadronic 
phase appears to come from the viscous Sf correction at 
freeze-out and thus not to be of dynamical origin. The 
sensitivity of the HBT variances to dynamical evolution 
effects caused by varying the temperature dependence 
of 77 /s at temperatures above freeze-out can be seen in 
Fig. 10 where the Sf contribution at freeze-out is re¬ 
moved. The results for viscous evolution split into two 
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FIG. 9: The normalized widths of the event-by-event distri¬ 
butions of HBT radii as functions of pair momentum Kt , 
for ideal fluid evolution and viscous dynamics with the four 
parametrizations shown in Fig. 1. 

bands, one for high, the other for low shear viscosity at 
early times. The relative widths for ideal fluid dynam¬ 
ics fall in between these two bands. The two bands are 
characterized by the behavior of 77 /s at high temperature, 
not at low T: Higher shear viscosity at early times (i.e. 
at high temperatures) leads to HBT variances that are 
10 - 20 % larger than those for evolution with lower shear 
viscosity at early times. The behavior of the shear vis¬ 
cosity at late times in the hadronic phase has no visible 
dynamical effect on the normalized widths of the HBT 
radii. 


IV. CONCLUSIONS 

We have presented a first analysis of the sensitivity 
of the HBT radii to the temperature dependence of the 
specific shear viscosity in relativistic heavy ion collisions 
at top RHIC energies, taking into account that the HBT 
radii fluctuate from event to event and that a consis¬ 
tent comparison of experimental data with theoretical 
predictions should thus be based on event-by-event evo¬ 
lution of fluctuating initial conditions on the theory side. 
We explored the effect of event-by-event fluctuations on 
the mean HBT radii (and their deviation from the val¬ 
ues obtained from single-shot hydrodynamics where the 
initial-state fluctuations are ignored by averaging over 
them before hydrodynamic evolution) and their relative 
widths. We found that fluctuations tend to increase the 
mean value of the fluctuating HBT radii by a few percent 
above the value obtained by the traditional process of av¬ 
eraging over the fluctuations already in the initial state, 
and found that the event-by-event distributions of the 
HBT radii are characterized by relative widths of order 
10 %. 


FIG. 10: Same as Fig. 9, but without the viscous correction 
Sf to the freeze-out distribution function. One sees that for 
the normalized width of the distribution of sideward radii R 2 , 
deleting Sf removes all sensitivity to the temperature depen¬ 
dence of 77 /s in the hadronic phase, while all other normalized 
widths are only weakly affected by Sf. 

Comparing the HBT radii from hydrodynamic evolu¬ 
tions with five different assumptions for the specific shear 
viscosity and its temperature dependence we showed that 
both the mean squared HBT radii and their relative 
widths are affected by a possible increase of the specific 
shear viscosity by a factor 10 between T tr and 3 T tr at 
the level of 10 - 20 %, with larger shear viscosities leading 
to smaller mean HBT radii and larger relative widths of 
their event-by-event distributions. There is little to no ef¬ 
fect on either from variations in the behavior of ( 77 /s) (T) 
in the hadronic phase below T tr , except for the sideward 
and out-side radius parameters R 2 and R 2 S 2 and the 
sideward variance a 2 which all show significant sensitiv¬ 
ity to the viscous Sf correction at freeze-out. The latter, 
in turn, depends on the value of 77 /s at the freeze-out 
temperature but not on its prior dynamical evolution in 
the hadronic phase. Since the Sf correction is not the¬ 
oretically well constrained at the present moment, this 
sensitivity of R 2 S and R 2 S 2 to 77 /s at freeze-out should be 
viewed with some caution. 

The levels of sensitivity of the HBT radii and their 
variances to the temperature dependence of 77 /s observed 
here are comparable to those observed in Ref. [1] for the 
elliptic and quadrangular flows. This indicates that pre¬ 
cise 3-dimensional HBT measurements, in particular new 
measurements of their variances [9], can play a valuable 
supporting role in constraining the temperature depen¬ 
dence of the QGP shear viscosity. 
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